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I. INTRODUCTION 

The main question we would like to study in this ar- 
ticle is how far can the notion of a particle, as a point 
like object be extended, in the realm of classical electro- 
magnetism. Of course in the textbooks one encounters 
extensive discussions of test particles; which are normally 
used for the very definition of the electromagnetic fields; 
namely, the electromagnetic fields are those that enter 
into the Lorentz force acting on test particles. However, 
the physical system becomes more complicated when one 
considers finite charged particles. On the one hand, such 
a particle radiates, and therefore the equations of motion 
should reflect the loss of energy-momentum of the parti- 
cle. On the other hand, a finite point like source involves 
fields that have a divergent behavior as one approaches 
the particle; which in turn would imply a divergent con- 
tribution to the total stress energy-momentum tensor of 
the system. 

This problem has been tackled in the past using dif- 
ferent techniques, as for example the idea of studying 
'spheres' in the limit when they become smaller and 
smaller. In fact this guided the early works of Lorentz, 
Abraham, and Dirac[T]. In this article we present a dis- 
cussion of the problem in which we avoid dealing with 
infinite terms. Our final result is a generalization of 
the equations found in the past, which contain previous 
works as particular choices. 

In the next section we present the basic notation used 
below. Section |III| is devoted to the presentation of the 
balance equations; which is our main tool in this work. 
The deduced equations of motion are a generalization of 
other similar equations found in the literature. We also 
present in this section a generalized notion of total mo- 
mentum for the charged particle. Notably the new equa- 
tions of motion involves new degrees of freedom, that are 
discussed below. A particular choice of the new degrees of 
freedom conduces to the celebrated Lorentz-Dirac equa- 
tions of motion. In section IV we study the nature of 
the solutions to the general problem. Those solutions 
present difficulties in the physical interpretation, which 



motivate us to study the equations of motion in terms of 
orders of the strength of the interaction with the external 
fields; which is presented in section |V] In the last section 
we summarize our results. An appendix is added where 
some properties of the solutions are presented. 



II. DEALING WITH PARTICLES IN 
ELECTROMAGNETISM 

A. Basic equations of a test charged particle 

Let us start by quickly reviewing the main equations of 
electromagnetism. The dynamics of the electromagnetic 
fields is governed by Maxwell equations; 



and 



(1) 



(2) 



While the dynamic of test particles is determined by the 
Lorentz force /" given by 



(3) 



where one is considering a test particle of charge q and 
four- velocity 

The immediate question is: what is the force that acts 
on a particle? (not test particle) 

The nature of the problem is that a (nontest) parti- 
cle radiates if it is accelerated; and therefore the use of 
the Lorentz force will imply an imbalance of energy and 
momentum. Several approaches have appeared to an- 
swer this question, and Lorentz, Abraham and Dirac[T] 
have made important contributions to this goal, which 
have ended in what is widely known as the Lorentz-Dirac 
equations of motion(see below). We here present these 
well known equations of motion as a special case of our 
equations by demanding balance of energy and momen- 
tum at an asymptotic sphere defined by the intersection 
of the future null cone of the particle with future null 
infinity. 
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B. Notation associated with a timelike curve 

Let us start with a massive point particle with mass 
771^ whose trajectory, in a flat space-time (M, 77^;,), is 
given by the timehke curve C, which in a Cartesian co- 
ordinate system reads 



(4) 



with r meaning the proper time of the particle along C . 
The associated 4- velocity of this particle is 



(5) 



and the signature of the flat metric is chosen such that 
VaV"" = 1. (Note that we are using units in which c = 1, 
so that w° has no units.) Now, for each point Q = z{t) 
of C, we draw a future null cone "^q with vertex in Q. If 
we call a:p the Minkowskian coordinates of an arbitrary 
point on the cone ^g, then we can deflne the retarded 
radial distance on the null cone by 



r^Va (a:?,-z"(T)). 
A null vector is defined by 

x° -z°(r) 
r 

And since Z° is null, one has 
Then one can see that 

vj'' = 1. 



(6) 



(7) 



(8) 



(9) 



One can introduce null polar coordinates associated 
to the timelike curve C in the following way. Let u be 
the null coordinate which is constant on the future null 
cones defined for points Q of C and such that u = t{Q). 
Let (C, C) be stereographic coordinates of the sphere of 
directions at the point Q. Then, the relation between 
Minkowskian coordinates and null polar coordinates 
(m, C: is given by 



= z''{u)+rr{uX.O- 



(10) 



Let us note that by defining Z° to be the null vector 
that corresponds to a stationary motion; which implies 
that P = ^°(C,C)j 0116 can see that Z° and 1°" wiU be 
proportional. In particular defining 



V{uX,O^V-{u)la{CX) 

one can see that 



UCX)- 



(11) 



(12) 



C. Electromagnetic fields of a nontest charged 
particle 

The retarded electromagnetic field of a particle with 
charge e, can be given in terms of the potential and/or 
the electromagnetic field itself. 

The retarded potential is 



ev°'(u) 



(13) 



The corresponding electromagnetic field is [2] 



ah 



2e 



2e 



V 



rV 



r^[aV6] 



(14) 



where a dot means derivative with respect to u and we 
have chosen Gaussian units; so that n, = An. 

It is convenient to also have at hand the spinor com- 
ponents of the electromagnetic field; which are given by 



0, 
e 

'2^' 



and 




(15) 



(16) 



(17) 



The symbol 3 denotes the edth operator^ of the unit 
sphere, and 9 is its complex conjugate. 



III. BALANCE EQUATIONS FOR A NON-TEST 
CHARGED PARTICLE 

A. The conservation law for the total energy 
momentum tensor 

The main tool to derive the balanced equations of mo- 
tion for charged particles is the conservation law of the 
total energy-momentum tensor; namely 



v-r = o. 



(18) 



Let be four translational Killing vectors; where we 
are using the numeric index c = 0, 1, 2, 3 and the abstract 
index b. Then from the conservation law ( 18 ) one also 
obtains that 



(19) 



Then, let y be the four volume which has as bound- 
aries the two spacelike hypersurfaces, E' at its future 



FIG. 1. Two spacelike hypersurfaces, E' and E, are the 
boundary of the internal four volume 'V . In turn, each of 
these spacelike hypersurfaces has the same two dimensional 
boundary S. 



boundary and E at its past boundary, as depicted in fig- 
ure [l] Also, each spacelike hypersurface has the same 
two dimensional boundary S. 

Given a three form D, Stoke's theorem tells us that 



dD 



D 



D. 



(20) 



In our case we take, for each Killing vector K^, the 
three form Dcabc = T'^e K^eabcd] so that dDcabcd = 
fcV/ (^T^^ K^^ £abcd', where A: is a constant. 
Therefore our main equation is 



= / dDc = 



Dr- 



(21) 



since due to the conservation equation ( 19 ) the left hand 
side vanishes. It is because of this reason that actually 
the total momentum is determined by S, and not by the 
particular hypersurface with boundary S. 

This equation is intimately related to the equations 
of motion; since equations (21) are telling us that the 



difference of the total momentum calculated at S' and 
S vanishes. We will apply these equations to the case of 
a particle, even if one has indications that each term in 
the difference could be ill defined; since it could contain 
infinite terms; however its difference is finite. 

Let iS be a sphere at future null infinity defined as 
the asymptotic sphere of the future null cone of a point 
Q{t); and let E be the future null cone of this point. Let 
S be the corresponding asymptotic sphere for the point 
Q{t + dr) on the curve C; and let E be the future null 
cone of this point. Let us call E+ the hypersurface at 
future null infinity with boundaries S and S. Then in 
the last equation we can identify E' = E U E+. See Fig. 
[2] for a graphical representation of the hypersurfaces. 

Therefore, assuming the same integrand, one has the 
relation 



(22) 



FIG. 2. Two hypersurfaces reaching future null infinity. 



which says that the momentum at S is the momentum 
at 5*, minus the flux through E+. 

Let us differentiate between a charged particle, which 
we will call the system of interest A, and the rest of the 
world which we will call system B. Then, the total energy 
momentum tensor can be decompose as the following sum 
of terms 



(m)A +r'(m)S 



-T, 



(EM)/1 " 



"T'(EM)B +T'(EM)A.b; (23) 



where we have distinguished: the mechanical term (m) of 
particle A, the mechanical term of system B, the electro- 
magnetic term (EM) of particle A, the electromagnetic 
term of system B, the electromagnetic term of the prod- 
ucts of fields of particle A times fields of system B. 

Considering the vector T^^K^; whose divergence is zero, 
one can apply equation (l22|). Each term of (23) con- 



tributes on the left hand side with a corresponding dP 
term, i.e. a difference of momentum at both null hyper- 
surfaces. The right hand side can be expressed as 



= - dr; 
s+ Js 



(24) 



where S is the intersection of the null cone S with future 
null infinity, and we are considering an infinitesimal dr. 
Therefore one has 



dP{ni)A + C^-P(m)B + dP(^EM)A + dP(^EM)B + 



dP, 



{EM)A. 



dP(EM)B.A^- J TabK^n-^dT 



dS^ , (25) 



where is a translational Killing vector, is the null 
vector normal to future null infinity satisfying n°'la = 1, 
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and dS'^ is the surface element of the unit sphere. We 
have also indicated separately the difference in momen- 
tum of system A due to the existence of system B 
(rfPfEMjA.s); E^nd the difference in momentum of system 
B due to the existence of system A (d-P(EM)s.A)- 

For the sake of simplicity in the discussion, let us con- 
sider first the case en which system B consists of a smooth 
distribution of matter and charges with no radiation at 
future null infinity. Furthermore, let us assume that sys- 
tem B can be represented by a Lagrangian formulation, 
included the actions of system A on B. Therefore, under 
variations of the B fields, the Lagrangian will induce the 
equations of motion 



dP, 



{m)B 



dPi 



(EM)S 



dP, 



dr 



dr 



in other words, the terms dPi 



(m)B - 



(EM)B.A 
-dPi 



0; (26) 



(EM)B " 

are balanced in equation (|25|), which implies 



dP, 



dP, 



(m)A 



dPi 



{EM)A 



dP, 



(EM)A,B 



dr 



dr 



dr 



{EM)B.A 



dS^ 



(27) 

Then one can think that the term (iP(EM)A could con- 
tain infinite contributions at the particle world line C, 
since the self fields will be taken into account. However 
we know that the right hand side of the equation is finite; 
therefore, the sum of both terms on the left is finite. 

To calculate the fiux term it is convenient to note that 



1 

47r 



dS^ 
y2 



(28) 



(29) 



and 



Now, let us note that for any two vectors A"' and B^, 
and defining the quantities A = A'^la and B = B°'la, one 
has 

A^Ba = AB+ABvSvB+BBvSvA~dvABvB~BvAdvB; 

(33) 

where we are assuming the Minkowski metric in the con- 
traction of vectors. 

Therefore, noting that y = v°'la, one deduces that 





where we have used that 



V 

V 



V 

v' 



(34) 



(35) 



Then, equation ( 32 1 becomes 




e 
47r 

2 



3 - 



(36) 



1 /■ ^92 1 

— / nH''— = 2w°i;^^ - - (rf^v'^ + if'v'' + r^'^'^v^) ; 
47r y 3 ^ ' ^ ' 

(30) 

which are a generalization of the relations that appear in 



The radiative part of Tab is given by 



1 



tab = 0202^a4; 



(31) 



Finally, the flux term contributes with 



TabKln'^ 



dS^ 



^ 2 -a • 

-e V VaVr. 

3 - 



(37) 



The other term on the right hand side of (271 includes 



all the contributions due to forces that system B exerts 
on particle yl; for the case of electromagnetic interactions 
one has 



so that the integration on S is more precisely 



2tt 



V \v -V 



(32) 



where 3y is the edth[3] operator for the metric with sur- 
face element 



dP, 



{EU)A,Bc 



= -eF{B)acv'' = eF{B),av''; (38) 



which is the standard Lorentz force. In order to elucidate 
the role of units, we clarify that we are using units of time 
for r (the proper time); so that a dot, or proper time 
derivative, implies an extra 1/second in any expression. 



The other terms on the left hand side of (27), not only 
include the mechanical term proportional to the acceler- 
ation Vc, but they also include other terms coming from 
the contributions due to the self fields of the particle. 
The main idea in our presentation is not to treat both 
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terms separately but as a unity; namely, the variation 
of the total momentum. In the customary treatment in 
which the terms are treated separately, one has to deal 
with infinities and complicated arguments for their can- 
cellation. We show here that the variation of the total 
momentum does not include difficulties. 

The variation of the total momentum does not have 
arbitrary vectorial dependence; since to begin with, the 
mechanical contribution only depends on the four veloc- 
ity, and the energy momentum tensor depends on the 
electromagnetic tensor, and this depends on the velocity 
and acceleration. Therefore, the time derivative of the 
total momentum can only depend at most on (v'^, i)'^, v'^); 
in other words it can not depend on higher derivatives 
like {'v'^). 

The coefficients of this vectors must be finite, since the 
right hand side is finite. Therefore, although at first sight 
one could expect divergent behavior for the expressions 
of the self fields; one must understand that there are 
cancellation of infinite like terms between the mechanical 
and the field momenta, that as a result provide with the 
finite coefficients. 

In other words, we do not attempt to decompose the 
variation of total momentum into a mechanical and elec- 
tromagnetic part, but consider an expression for the vari- 
ation of the total momentum; which we know is finite, 
and of the form 

Pc — m'vc + 'mvc + a:Vc, (39) 
were m', m and a are understood as functions of r. 



From equation (271 we conclude then that the equa- 



tions of motion can be written as 

Let us also note that the total momentum must be 
expressed as 



Pr = MVr 



avc 



(41) 



with M = M{t). This expression for the total momen- 
tum is a generalization of the expression (4.4) found in 
reference [S] through renormalization procedures. 

By taking the time derivative of the last equation and 
comparing with eq. ( 39 ) , we get 



M = to', 
M + a = m. 



(42) 
(43) 



Let us note that there are at least two notions of mass. 
One has the parameters to, which we will call inertial 
mass. Also one has the parameter M, that we will call 
rest mass, in order to differentiate them. 

Therefore the final form of the equations of motion is 



Mvc + (M + a)vc + aVc = eF{B)caV°- + -e^v^VaVc 



This is our main result, the most general equations of 
motion for charged particles in the framework we have 
presented; in which it is required the balance of radi- 
ated momentum. Let us note that in our approach we 
have avoided treating explicitly the elimination of infinite 
renormalizations or similar techniques. 

We recognize new degrees of freedom present in this 
form of the equations of motion, namely the scalar Q;(r); 
and also M(r). The latter has been consider already by 
other authors (see below). 

Now, if we contract the expression eq.(44) with the 
vector and using the notation — = v°'Va, we obtain 



M 



' 2 2 

-e a . 



(45) 



Since this is one equation for two unknown, one must 
prescribe o.{t) or M(r), or a relation among them, as we 
will do below. 

It is interesting to note that if we replace M in the last 
equation into ( 44 ) one obtains 



(A/ + a)vc = eF{B), 



(46) 



which only sets the dynamics for three degrees of freedom 
since it is orthogonal to the four velocity v. Therefore one 
has the option to consider the original equations ( 44 1 , or 
the two independent equations ( 46 1 and ( |45l). 

In presenting equations ( 44 1 , ( 45 1 and ( 46 1 we have 



given preference to the rest mass M; alternatively, if one 
gives preference to the inertial mass m one can reexpress 
them as 



(TO - aw + mv, + av, = eFiBUv'^ + -e^v'^v^v. 



m — a + ac? = — e^a^, 
3 



and 



(47) 
(48) 

(49) 



(44) 



In summary, we have arrived at the general equations 
of motion for a charged particle, in which two new degrees 
of freedom appear, and also the order of the equations 
for V has changed; from the Lorentz equations. So, the 
set of basic variables, can be considered {v°',M,a) or 
(ti", TO, a); where the velocity vector is in turn expressed 
in terms of the position of the particle. Since v is always 
assumed to have unit modulus, one has four equations for 
five degrees of freedom. Therefore one has the liberty to 
choose a relation among the two new degrees of freedom. 
We will discuss below how different choices conduces us 
to known cases and also to new physically interesting 
ones. 

It is probably worthwhile to emphasize that we have 
arrived at these general equations of motion by balanc- 
ing the retarded radiation field of the charged particle 
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when the rest of the system does not radiate. This is a 
very strong assumption, on the rest of the electromag- 
netic system, and in particular would exclude the inter- 
action of the particle with perfect conductors, or a dis- 
persive permeable medium, for example. Our approach 
is a technique valuable to calculate the 'corrections' to 
the Lorentz force that take into account the retarded 
radiation fields; however, after one has calculated such 
corrections one expects these to be valid in a general sit- 
uation. 



Historical choices on the general equations of 
motion 



1. The Lorentz-Dirac equations 

We have seen before that there are two natural notions 
of mass that appear in our treatment. Therefore it is 
interesting to consider the particular cases in which each 
of them is required to be a constant of the motion; since in 
particular real elementary particles as the electron seem 
to have a constant mass. 

The case m = will be considered below; we will con- 
sider here the case M = 0. 

This particular choice can be understood as follows. 
From (411, we see that M — m — a can be interpreted 



as the rest mass of the electron (which has contributions 
of the electromagnetic field generated by the particle). 
Then if we want to describe particles whose total rest 
mass remains constant, we must require that M = rh — 



d = 0, then from eq.(45l, we obtain that 



(50) 



Furthermore, one can also deduce from equation ( 43 ) 



that the two notions of mass coincide, and therefore one 
also has rh = 0. 

Therefore, one arrives at the following equations of mo- 
tion for the charge given by 



mv 



eF{Bry + -e^r + v''vhv''); (51) 



which are the well known Lorentz-Dirac equations of mo- 
tion for the electron. 

It is probably worth while to recall that these equa- 
tions have several problems, as it has been investigated 
in the past by several authors. First of all they are third 
order differential equations for the position of the charged 
particle; which is contrary to the general accepted idea of 
mechanics for particle. Secondly, they have the so called 
problem of the runaway solutions. Several of the implied 
problems were discussed by Dirac and we will not review 
them here. 



2. The Bonnor equations 

Another approach to the two mass question is to re- 
quire that the two notions coincide. Then, let us consider 
the choice M = m. From this one immediately obtains 
that d = 0; and we recognize that the choice a = — 
coincides with the previous case; but instead we study 
here the general rh ^ case; which must satisfy 



m 



(52) 



where as noted before, a must be a constant. 

In particular if we set a = one arrives at the equa- 
tions of motion with varying mass, namely 



raVc 



with 



TO 



' 2 2 

-e a . 



(53) 



(54) 



This equation was studied by |6j. Unfortunately it has 
the unphysical consequence that the mass of the particle 
could vanish in a finite time. 



IV. STUDY OF THE GENERAL EQUATIONS 
OF MOTION FOR CHARGED PARTICLES 

A. Behavior of the general solution for two 
particular cases 

We study here general properties of the solutions that 
can be deduced prior to the use of the extra liberty for 
the choice of a condition among the two new degrees of 
freedom. 

a. The case a" — Q. Let us consider first the case in 
which a? = for all times. In this case, all the terms in 



equations (44) collapses to zero, and the particle moves 



along a geodesic. One can observe that the terms involv- 
ing i) and v are zero; then the only two remaining terms 
containing the external fields and the first one propor- 
tional to V are orthogonal, so that they must vanish in- 
dependently. In particular one has F{B) = and M = 0. 

h. The case F{B) = 0. The other case that it is im- 
portant to be considered is the case in which the external 
fields F{B) are zero, and study the general behavior of 

for large values of the time variable; in order to see the 
nature of the runaway solutions problem in this setting. 
So now we invert the logic and think what is the behav- 
ior of in terms of the behavior of the new degrees of 
freedom M and a. We observe then that by contracting 
the equations of motion ( 44 ) with v one obtains 



{M + a)a^ + |a2 = 0. 



(55) 



This constitutes a simple first order differential equation 
for a^; whose behavior is determined by the sign of the 
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coefficients. In this setting a runaway behavior imphes 
that the coefficients have opposite sign. Since one can 
choose one relation for the two extra degrees of freedom, 
we would like to explore possible conditions for them. If 
one were interested in an asymptotic nonincreasing solu- 
tion for a^, one would probably consider {M + a) > 
and a > 0. It is interesting to note that if one chooses 
as initial data for a and a a tiny positive value for them. 



then equation ( 55 1 would imply an initial fast exponen- 
tial decay of a . This , will in turn imply, that initially 
M ^ m, and from (45) that M would become very small 
in a short time. This choice of very small initial condi- 
tions of the a degree of freedom avoid the problem of the 
runaway solutions, for the case in which F{B) = for all 
times. Instead, the discussion becomes more complicated 
if one thinks in the situation in which the external field is 
turned on and off. When the external field is turned on, 
one can arrange so that > 0; which will imply from 
(551 that a rapidly adopt a negative value (provided one 
has manages to maintain m — M + a > 0). If later the 
external field is turned off again, then the particle will 
enter the region of vanishing field with the wrong initial 
condition and one would probably be in the presence of 
a runaway solution. So there is no universal choice of the 
initial conditions of the a degree of freedom that would 
exclude the runawa y so lutions. Also, the behavior a > 
would imply, from (45), a condition M < 0, which could 



have the problem that M would vanish in a finite time, if 
where bounded from below by a nonzero value. From 
this analysis we can not exclude the possibility of a solu- 
tion where both a and go asymptotically to zero, for 
large values of the time coordinate, but with a nonzero M 
value. The other possibility is to have a < 0; but would 
conduce us to the runaway solution problem (provided 
one has manages to maintain m = M + a > 0). There- 
fore, we conclude that the runaway solution problem is 
generic in this dynamics too. The case (M -I- d) < is 
ruled out since this would imply that the inertial mass is 
negative. The value a = was not considered due to the 
fact that it would collapse to the previously studied case 
of = (since we are examining the case F{B) = 0). 



B. Behavior of the solution with m = 

Let us recall that m can be understood as the iner- 
tial mass. The reason that condition m = is worth 
studying comes also from the following considerations. 
The time derivative of the rest mas^ of the particle is 
given by M = TO — d. Therefore, the condition m — 
means that the inertial mass is constant, and the varia- 
tion in the rest mass is only due to the electromagnetic 
fields contribution to the mass. Note that the Lorentz- 



^ Recall that M is the factor of v in the expression for the total 
momentum. 



Dirac equations are contained in this family; as it is case 
F{B) = considered in the last subsection. 

Treating the equations of motion as exact to any or- 
der leads one to possible divergent behavior for the new 
degrees of freedom as we will see next. 



In this case, equation (45 1 can be completely expressed 



in terms of a; which must satisfy 



2 2 2 

a — a a = -e a . 



(56) 



It is convenient to define /3 = a+ |e^, since then eq.(56l, 
can be written as 



/3- 



'/3 = 0. 



(57) 



In appendix [A] there is a brief discussion of the prop- 
erties of solutions of equation ( 57 ) in terms of the global 
properties of a^(T). There, a list of theorems is pre- 
sented where nondivergent properties for are assumed 
for large values of r. These results suggest that one could 
in principle be able to choose initial conditions for /3 to 
select decaying solutions. That is, one in general would 
find independent solutions for /3; one that would grow 
with time and other that would decay to zero. It is some- 
how remarkable that one could select a solution for f3 so 
that a tends asymptotically to the Lorentz-Dirac value. 
However, we emphasize that in order to find the decreas- 
ing solution one must know the whole future history of 
the motion of the particle; therefore this is a nonlocal 
analysis of the motion; which in particular is not based 
just on initial conditions. 

Let us note that defining the vector — eF{BY^^v'', 
equations ( 49 1 become 



mv 



(58) 



which with (56) constitute two independent equations. 



Equation (58) has the same form as the Lorentz-Dirac 



equations; and it would agree with it if one would take 



the solution a = 



2„2 



as mentioned above. 



There is a general concern with the dynamical equa- 
tions chich involves the choice of two constants for the 
initial conditions for equation (56). Then in turn one 



would have the set of solutions of the equations of mo- 
tion ( 44 ) affected by the arbitrary initial choice for the a 



degree of freedom. 



V. DYNAMICS IN TERMS OF ORDERS OF 
THE STRENGTH OF THE INTERACTION 
WITH EXTERNAL FIELDS 

The difficulties found in the study of the general solu- 
tions, of the last section, suggest that the original equa- 
tions of motion must be understood in terms of an anal- 
ysis in orders of the interaction with the external fields, 
which we present next. 
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Equation of motion up to second order in the 
external fields 



It is important to remark that equations ( |44[ ) are exact 
equations in the classical framework for charged parti- 
cles. It is generally believed that a classical description 
of particles must break down when one considers real ele- 
mentary particles at microscopic levels. Presumably, one 
would be able to describe the behavior of the new degrees 
of freedom from quantum electrodynamics. All this sug- 
gests that equations ( 44 1 should be understood in terms 



of orders of the interaction with the external field F{B). 
Therefore next, we present a study of our main equa- 



tions ( 44 ) in first orders of the strength of the interaction 



with the external field; but still considering the physical 
condition rh — 0. This condition allows to interpret M(t) 
in terms of a^r). 

Note that contracting /° with eq. (44), one obtains 



mrva + arva = -r, 



(59) 



with f'^ = — /q/". At this point it is convenient to re- 
call that the equations of motion for a test particle with 
negligible mass and charge /i and q, respectively, is 



^J.Vc 



qF{B),av''; 



(60) 



from which one deduces that the behavior of a{T) is com- 
pletely due to the interaction of the charge e with the 
external field Fab- If one considers the correction to the 
equations of motion as arising from the strength of dif- 
ferent terms appearing in eq. ( 44 ) , one is tempted to con- 
sider terms of order ^(P) and also to consider order in 
^{e'^); since the appears as an independent factor in 
the radiation term. 



Let us note that from equations ( 60 1 one can deduce 
that V = GiX)- Then, since we know that a - 



would have, from equation (59 1 that v — — ^(f) 



Then, since a is expected to have some nontrivial i^(e'') 



order, the second term on the left hand side of ( 59 1 is of 



higher order than the rest; as can also be deduced from 
the fact that 



1 



m 



(61) 



where ^(f"*") means higher order than (f), as for example 
^(e^f), and for simplicity we have omitted the vectorial 
index. 

Now we will consider equations (58) up to order ^(f^). 
We therefore study equation (57) up to order ^^(f^). For 
this purpose we will use a and suggest a /3 of the 

form 

(3{t) = Ao(t) + Ai(r)f + A2(T)f2 + ^(f3). (62) 



Then equating order like terms in equation (57) one 
has: 



Ao = 0, 



(63) 



and 



(Aif) = 0, 



(A2f2) = Ao 



(64) 



(65) 



The solution of equation ( 64 ) implies two constants of 



integration that must be of ^(f); but for the general 
physical situation of a nonstationary f we find no univer- 
sal way to assign a constant. Therefore we chose those 
two constants to be zero; in other words, we set Aif = 0. 



Instead in the solution of (63) one finds Aa = Aqo 
Aqit; where Aqq and Aqi are constants. 



Replacing the solution of these terns in equations ( 58 ) 
and keeping terms up to if5'(f^) one obtains 



mv'^ = r + i^e' - Ao) ( -r 



(66) 



It is hard to give a physical meaning to a dynamical sys- 
tem that depend on the two parameters involved in Aq. 
Even if one chooses Aqi = to avoid the time depen- 
dence; one is still left with a constant that probably 
should be determined from another theoretical frame- 
work like quantum electrodynamics. From the classical 
point of view one could consider studying the cases in 
which Aq is taken to be a constant proportional to . A 
peculiar case would be to consider Aq = |e^; since this 
choice would just cancel the other terms proportional to 
e^, and therefore conduce to the Lorentz force for the 
charged particle. This would be unacceptable in our ap- 
proach to the equations of motion that takes into account 
the change in momentum due to the radiation emitted by 
the particle. However, considering a background quan- 
tum nature of particles, one could think in a value of the 
form 



Ao^X- 



(67) 



where h is Planck constant, c the velocity of light and A 
a number without units that should be determined from 
quantum electrodynamics. Then, since we are consider- 
ing second order expressions, in the equations of motion 
( 66 ) , in place of fc one must use 



/(2)c ={eF{B),t. yb) = e(i^(B),fc) v'' + eF{B),h 

=e{F{B),b) v'>+-F{BUf. 

m ^ 

(68) 

In a complete classical framework, one would not have 
any physical argument for a nonzero Aq\ and therefore 
one would be forced to take Aq = 0; leading to the equa- 
tions of motion 



6 m ^ ' 6 



(69) 



which coincides with the Landau-Lifshitz[7j equations of 
motion. 
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B. Study of the equations of motion up to third 
order in the interaction with external fields 



Now we would like to study equations ( 44 1 , or equiva- 



lently (49) up to third order. 



Let us note that by having obtained the second order 
equations of motion in the previous section; one has a 
second order acceleration ^(2); namely, the one that sat- 
isfies equations ( 69 ) . Then, third order in turn is defined 



to satisfy the equations of motion 



mv 



(3) 



"(2)(W('2) 



^(2) 



(70) 



where ii'l^^) ^^e r derivative of 

«^2) = -f + ile'-A,)^J^,^-{le'-Ao)^^v'; (71) 

and one should consider terms up to order 3. 

In this way one can see that this procedure can be 
generalized to any higher order where satisfies the 
differential equations 



All this seems to indicate that the only physically sen- 
sible choice is to take Aq = 0; since in this case one 
would have to solve the homogeneous problem, with the 
natural choice for all arbitrary constants to be zero. In 
other words, /3 must be zero and we are conduced to the 
equations 



'(3) 



^(^(2) 



where 



2 o 1 
A m 



2 /(2) 



2 . 

ttG — 5-u ; 



(77) 



(78) 



in which we have replaced a(2) = ^f^^ one should 
consider terms up to order 3. 

This difficulty will appear in higher orders too; from 
which we conclude that a(„^2) — "§6^, and that the 
second order case treated previously is sort of peculiar, 
since it is the only order that allows for a physically per- 
missible nonzero /?. 



) rb 



(72) 



VI. FINAL COMMENTS 



where ii/ ,n is the time derivative of 



V(n-1) = - "(«-2) (^('n-2) " a(„_2)^^''); (73) 

and one should only consider terms up to order n. 

Then, coming back to the third order calculation, let 
us consider /3 of the form 

/3(r) = Ao(T)-|-Ai(r)f+^2(r)f2+A3(r)f3 + ^(f4). (74) 



Then equating order like terms in equation (57) one ob- 



tains equations (63) and (64); while instead of (65) one 
now has 



and also 



(A3f3) 



Aoo^ 



(3)' 



(75) 



(76) 



where a?^^) and a^^., are the order i^(f^) and ^(f^) re- 



"(3) 



spectively of . 

If ^ 0, the main difhculty in this case is the fact 
that in the product Pa? of equation (57) one needs the 



explicit integral of (75), which will involve a couple of 
integration constants(or order i^(f^)) associated to the 
choice of initial time for the integration. Since the inte- 
gration constants must be of order ^(f^), they can not 
be associated just with the charge 'e' or to universal con- 
stants; therefore it would be very difficult to give physical 
meaning to a physical dynamical system that depends on 
the arbitrary choice of initial conditions (even if one man- 
age to choose them of order ^(f^)). 



The equations of motion ( 44 ) are general equations for 



charged particles that are derived from the condition of 
balance of variation of total momentum with the radiated 
momentum. In this derivation, we have avoided dealing 
explicitly with infinite contributions to the momentum. 

We have shown how different versions of equations of 
motion for charged particles can be obtained from our 
general equations; in particular the celebrated Lorcntz- 
Dirac equations. 

Studying the properties of solutions to the general 
equations of motion we have found the possibility to 
choose initial conditions for the new degrees of freedom, 
which select decaying modes to the Lorentz-Dirac value. 
However this choice is possible only if one knows the 
whole history of the world line of the particle. This is 
reminiscent of the notion of a horizon in general relativ- 
ity, which it depends on the whole history of the space- 
time. Even if one manage to chose these preferable initial 
conditions, one is still left with the problem of runaway 
solutions, or equivalently of pre-accelerations. If one were 
forced to a theoretical framework based purely on clas- 
sical considerations, this would be a sort of dead-end. 
Nonetheless, if one uses the known physical information 
that the ultimate nature of real particles is of a quan- 
tum kind; and therefore classical theoretical frameworks 
should be understood as approximate models of the real 
world, one would think that one should not demand a real 
physical interpretation to the exact classical equations. 
Instead one is tempted to consider the corrections to the 
Lorentz force, as also being related to terms that proba- 
bly should be calculated from quantum electrodynamics. 
If so, then the original equations should be understood 
in terms of orders of the strength of the interactions with 
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the external fields. In particular there is a natural limit 
to the strength of the interactions that must be consid- 
ered in order to avoid the quantum creation of pairs of 
particles. 

In the study in terms of orders of the strength of in- 
teraction with the external fields we have found at sec- 
ond order the equations of motion ( 69 1 . These equations 



Then the equation 



have also been supported in the derivation of reference 
[S] based on a different setting; namely, the study of first 
effects for particles with 'small' charge and mass. Instead 
herewe have considered finite charge and mass particles 
but study the exact equations of motion in terms of or- 
ders of the strength of the interaction with the external 
fields. Then, it is no surprise that both approaches agree 
at first orders. 

It is important to remark however that the original 
Landau-Lifshitz equations are of second order and coin- 
cide with our equations ( 69 ) . Our second order equations 



( 66 1 could account for first order quantum corrections to 



the classical equations. 

Instead, the equations of motion (72) and (73) pre- 



sented here, are a generalization valid up to any desired 
higher order. 

Summarizing; although we have found more general 
equations of motion for charged particles, we have shown 
that our set contains the main cases studied in the past. 
We also give arguments that indicate that one should not 
take the exact general equations (44) as the physically 
relevant ones. Instead, to our understanding the only 
physically reasonable treatment of the equations of mo- 



tion ( 44 ) is through the notion of finite orders in terms of 



so that we conclude that the equations of motion 
applicable to classical particles, but with finite charge, 
are (|69]), in second order, or equations (72) and (73) in 
higher orders. 

The final appendix section presents properties of the 
solutions to equation (57). 
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Appendix A: Properties of /3 function 

Here we would like to recall some of the properties of 



the solutions to equation (57) 



Th. A.l|^ Let a^(T) > be continuous in (—00,00). 



(57) 



has one and only one solution /3i{t) passing through (0,1) 
which is positive and strictly decreasing for all r and one 
and only one solution /32(t) through (0,1) which is positive 
and strictly increasing for all t. Furthermore 

a^i e L2(0, 00), Pi e £2(0, 00). 

The notation for ^2(0, 00) comes from the Lebesgue 
spaces Lp{b, c); which contain the set of integrable func- 
tions such that 



< 00. 



(Al) 



Th. A. 2 If < a?f < < a^^jp < 00, then the solution 
(3i satisfies 



< /3(t) < e--^-^; 



(A2) 



for T > 0. 



Th. A. 3 I^Let be positive and continuous in [0, 00) and 
posses continuous first and second order derivatives. Set 



H = 



^(„2)2_ 1„2„2 



and suppose that 

H e Li(0,oo). 
Then there exist constants ci and C2 such that 

— pIo a{T')dT' 

/3i(r) = ci ^ (1 + i?i(T)) 

\/a 



(A3) 
(A4) 

(A5) 



efo 

/?2(r)=C2 ^ (l+i?2(T)), (A6) 



where 



|i?i(T)K eJ'-°°l^("')l'^"' - 1; (A7) 
and a positive constant c such that |i?2(T)| ^ c|i?i(T)|. 

Th. A.4 0ln the equation 

/3 - /? = 0, 



(57) 



let be a continuous complex-valued function for large t 
satisfying 



/oo 
r|a^(T)|dT < 00, 



(A8) 



^ Page 359 in reference [9]- 



^ Page 445 in reference [9]. 
* Page 380 in reference [10] 
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or, more generally 



{T')dT' = lim J a^{T')dT' exists and 

/oo 
sup |A(r)|(ir < oo. 
r<r<oo 



Then, there is a pair of solutions j3o and /3i satisfying, as 

T — )■ OO, 



/3o(t)^1, /3o(t)=o 



/3i(T)^r, /3i(r)^l. 



(A9) 



(AlO) 



Conversely, if a?{T) is real-valued and does not change signs 
and if ([57]) has a solution satisfying ( [A9| or ( |A10[ ), then ( [AS] ) 
holds. 



Th. A.5 0ln the equation 

/3 - (A2 + g(T)) /3 = 0, 



(All) 



let A > and q{T) be a complex-valued continuous function 
for large r satisfying 



|g(T)|dT < 00, 



(A12) 



or, more generally. 



lini 

T-i-oo 



e^^^ sup 



exists and 



g(r)e 



(A13) 



dr < 



Then, (All) has solutions /3o, /^i satisfying 



/3o 



' A 



(A14) 



Conversely, if (/(r) is real-valued and does not change signs 
and if (All) has a solution (3o or Pi satisfying the corre- 



sponding conditions in (A14|, then (A12| holds 



An example of a solution that does not satisfy any of 
the hypothesis of the theorems is the case = The 
two independent sohitions for /3 are proportional to 
and to l/r. 

However, let us note that the case 



satisfy the hypothesis of the last theorem; which states 
that there is an exponentially decreasing solution of the 
form /?o ^ e'^" . 
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